Search Results

Simultaneous two-component measurements of gas velocity and multi-dimensional numerical simulation are employed to characterize the evolution of the in-cylinder turbulent flow structure in a re-entrant bowl-in-piston engine under motored operation. The evolution of the mean flow field, turbulence energy, turbulent length scales, and the various terms contributing to the production of the turbulence energy are correlated and compared, with the objectives of clarifying the physical mechanisms and flow structures that dominate the turbulence production and of identifying the source of discrepancies between the measured and simulated turbulence fields. Additionally, the applicability of the linear turbulent stress modeling hypothesis employed in the k-ε model is assessed using the experimental mean flow gradients, turbulence energy, and length scales.

A level set G-equation model has been developed to model the combustion process in spark ignition engines. The spark ignition process was modeled using an improved version of the Discrete Particle Ignition Kernel (DPIK) model. The two models were implemented into the KIVA-3V code to simulate SI engine combustion under both premixed and direct injection conditions. In the ignition model, the ignition kernel growth is tracked by Lagrangian markers, and spark discharge energy and flow turbulence effects on the kernel growth are considered. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the G-equation model is used to track the mean turbulent flame evolution. When combined with a characteristic time scale combustion model, the models were also used to simulate stratified combustion in DISI engines, where the triple flame structure must be considered.

An improved discrete particle ignition kernel (DPIK) model and the G-equation combustion model have been developed and implemented in KIVA-3V. In the ignition model, the spark ignition kernel growth is tracked by Lagrangian markers and the spark discharge energy and flow turbulence effects on the ignition kernel growth are considered. The predicted ignition kernel size was compared with the available measurements and good agreement was obtained. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the level set method (G-equation) is used to track the mean turbulent flame propagation. It is shown that, by ignoring the detailed turbulent flame brush structure, fine numerical resolution is not needed, thus making the models suitable for use in multidimensional modeling of SI engine combustion.

The present study uses a numerical model to investigate the effects of flow turbulence on premixed iso-octane HCCI engine combustion. Different levels of in-cylinder turbulence are generated by using different piston geometries, namely a disc-shape versus a square-shape bowl. The numerical model is based on the KIVA code which is modified to use CHEMKIN as the chemistry solver. A detailed reaction mechanism is used to simulate the fuel chemistry. It is found that turbulence has significant effects on HCCI combustion. In the current engine setup, the main effect of turbulence is to affect the wall heat transfer, and hence to change the mixture temperature which, in turn, influences the ignition timing and combustion duration. The model also predicts that the combustion duration in the square bowl case is longer than that in the disc piston case which agrees with the measurements.

Simultaneous measurements of the radial and the tangential components of velocity are obtained in a high-speed, direct-injection diesel engine typical of automotive applications. Results are presented for engine operation with fuel injection, but without combustion, for three different swirl ratios and four injection pressures. With the mean and fluctuating velocities, the r-θ plane shear stress and the mean flow gradients are obtained. Longitudinal and transverse length scales are also estimated via Taylor's hypothesis. The flow is shown to be sufficiently homogeneous and stationary to obtain meaningful length scale estimates. Concurrently, the flow and injection processes are simulated with KIVA-3V employing a RNG k-ε turbulence model. The measured turbulent kinetic energy k, r-θ plane mean strain rates ( 〈Srθ〉, 〈Srr〉, and 〈Sθθ〉 ), deviatoric turbulent stresses , and the r-θ plane turbulence production terms are compared directly to the simulated results.

An optimization study combining multidimensional CFD modeling and a global, evolutionary search technique known as the Genetic Algorithm has been carried out. The subject of this study was a 2-stroke, spark-ignited, direct-injection, single-cylinder research engine (SCRE). The goal of the study was to optimize the part load operating parameters of the engine in order to achieve the lowest possible emissions, improved fuel economy, and reduced wall heat transfer. Parameters subject to permutation in this study were the start-of-injection (SOI) timing, injection duration, spark timing, fuel injection angle, dwell between injections, and the percentage of fuel mass in the first injection pulse. The study was comprised of three cases. All simulations were for a part load, intermediate-speed condition representing a transition operating regime between stratified charge and homogeneous charge operation.

A radiation model based on the Discrete Ordinates Method (DOM) was incorporated into the KIVA3v multidimensional code to study the effects of soot and radiation on diesel engine performance at high load. A thermophoretic soot deposition model was implemented to predict soot concentrations in the near-wall region, which was found to affect radiative heat flux levels. Realistic, non-uniform combustion chamber wall surface temperature distributions were predicted using a finite-element-based heat conduction model for the engine metal components that was coupled with KIVA3v in an iterative scheme. The more accurate combustion chamber wall temperatures enhanced the accuracy of both the radiation and soot deposition models as well as the convective heat transfer model. For a basline case, (1500 rev/min, 100% load) it was found that radiation can account for as much as 30% of the total wall heat loss and that soot deposition in each cycle is less than 3% of the total in-cylinder soot.

Optimizations were performed on a Heavy-Duty diesel engine equipped with a conventional electronic unit injector in order to minimize fuel consumption, and emissions of NOx and particulate matter. A low speed light load case and a high speed light load case were optimized with these considerations in mind. Exhaustive parametric studies were performed in order to find sets of operating conditions that resulted in low emissions and high fuel economy. It was found for the low speed light load case (Mode 2, 25% load and 821 rev/min) that low emissions operating conditions existed at either very early or very late start-of-injection timings and high EGR (PM = 0.018 g/kW-hr, NOx + HC = 1.493 g/kW-hr with SOI = -21 degrees ATDC, 48% EGR; or 0.085 g/kW-hr PM, 1.02 g/kW-hr NOx with SOI = 4 degrees ATDC, 39% EGR).

To overcome the trade-off between NOx and particulate emissions for future diesel vehicles and engines it is necessary to seek methods to lower pollutant emissions. The desired simultaneous improvement in fuel efficiency for future DI (Direct Injection) diesels is also a difficult challenge due to the combustion modifications that will be required to meet the exhaust emission mandates. This study demonstrates the emission reduction capability of split injections, EGR (Exhaust Gas Recirculation), and other parameters on a High Speed Direct Injection (HSDI) diesel engine equipped with a common rail injection system using an RSM (Response Surface Method) optimization method. The optimizations were conducted at 1757 rev/min, 45% load. Six factors were considered for the optimization, namely the EGR rate, SOI (Start of Injection), intake boost pressure, and injection pressure, the percentage of fuel in the first injection, and the dwell between injections.

The present study investigated HCCI combustion in a heavy-duty diesel engine both experimentally and numerically. The engine was equipped with a hollow-cone pressure-swirl injector using gasoline direct injection. Characteristics of HCCI combustion were obtained by very early injection with a heated intake charge. Experimental results showed an increase in NOx emission and a decrease in UHC as the injection timing was retarded. It was also found that optimization can be achieved by controlling the intake temperature together with the start-of-injection timing. The experiments were modeled by using an engine CFD code with detailed chemistry. The CHEMKIN code was implemented into KIVA-3V such that the chemistry and flow solutions were coupled. The model predicted ignition timing, cylinder pressure, and heat release rates reasonably well. The NOx emissions were found to increase as the injection timing was retarded, in agreement with experimental results.

Cycle-resolved analysis of velocity data obtained in the re-entrant bowl of a fired high-;speed, direct-injection diesel engine, demonstrates an unambiguous, approximately 100% increase in late-cycle turbulence levels over the levels measured during motored operation. Model predictions of the flow field, obtained employing RNG k-ε turbulence modeling in KIVA-3V, do not capture this increased turbulence. A combined experimental and computational approach is taken to identify the source of this turbulence. The results indicate that the dominant source of the increased turbulence is associated with the formation of an unstable distribution of mean angular momentum, characterized by a negative radial gradient. The importance of this source of flow turbulence has not previously been recognized for engine flows. The enhanced late-cycle turbulence is found to be very sensitive to the flow swirl level.

In-cylinder fluid velocity is measured in an optically accessible, fired HSDI engine at idle. The velocity field is also calculated, including the full induction stroke, using multi-dimensional fluid dynamics and combustion simulation models. A detailed comparison between the measured and calculated velocities is performed to validate the computed results and to gain a physical understanding of the flow evolution. Motored measurements are also presented, to clarify the effects of the fuel injection process and combustion on the velocity field evolution. The calculated mean in-cylinder angular momentum (swirl ratio) and mean flow structures prior to injection agree well with the measurements. Modification of the mean flow by fuel injection and combustion is also well captured.

A spray model for pressure-swirl atomizers that is based on a linearized instability analysis of liquid sheets has been combined with an ignition and combustion model for stratified charge spark ignition engines. The ignition model has been advanced, such that the presence of dual spark plugs can now be accounted for. Independent validation of the spray model is achieved by investigating a pressure-swirl injector inside a pressure bomb containing air at ambient temperature. In a second step, the complete model is used to estimate the performance of a small marine DISI Two-Stroke engine operating in stratified charge mode. Simulation results and experimental data are compared for several different injection timings and the agreement is generally good such that there is confidence in the predictive quality of the combustion model. Finally the model is applied in a conceptual study to investigate possible benefits of split injections.

Detailed chemical kinetics was implemented in the KIVA-3V multidimensional CFD code to study the combustion process in Homogeneous Charge Compression Ignition (HCCI) engines. The CHEMKIN code was implemented such that the chemistry and flow solutions were coupled. Detailed reaction mechanisms were used to simulate the fuel chemistry of ignition and combustion. Effects of turbulent mixing on the reaction rates were also considered. The model was validated using the experimental data from two modified heavy-duty diesel engines, including a Volvo engine and a Caterpillar engine operated at the HCCI mode. The results show that good levels of agreement were obtained using the present KIVA/CHEMKIN model for a wide range of engine conditions, including various fuels, injection systems, engine speeds, and EGR levels. Ignition timings were predicted well without the need to adjust any kinetic constants.

A new ignition and combustion model has been developed and tested for use in premixed spark-ignition engines. The ignition model is referred to as the Discrete Particle Ignition Kernel (DPIK) model, and it uses Lagrangian markers to track the flame-front growth. The model includes the effects of electrode heat transfer on the early flame kernel growth process, and it is used in conjunction with a characteristic-time-scale combustion model once the ignition kernel has grown to a size where the effects of turbulence on the flame must be considered. A new term which accounts for the effect of air-fuel ratio, was added to the combustion model for modeling combustion in very lean and very rich mixtures. The flame kernel size predicted by the DPIK model was compared with measurements of Maly and Vogel. Furthermore, predictions of the electrode heat transfer were compared with data of Kravchik and Heywood. In both comparisons the model predictions were in good agreement with the experiments.

A correction for the turbulence dissipation, based on non-equilibrium turbulence considerations from rapid distortion theory, has been derived and implemented in combination with the RNG k - ε model in a KIVA-based code. This model correction has been tested and compared with the standard RNG k - ε model for the compression and the combustion phase of two heavy duty DI diesel engines. The turbulence behavior in the compression phase shows clear improvements over the standard RNG k - ε model computations. In particular, the macro length scale is consistent with the corresponding time scale and with the turbulent kinetic energy over the entire compression phase. The combustion computations have been performed with the characteristic time combustion model. With this dissipation correction no additional adjustments of the turbulent characteristic time model constant were necessary in order to match experimental cylinder pressures and heat release rates of the two engines.

A study was performed that takes into account both turbulence and chemical kinetic effects in the numerical simulation of diesel engine combustion in order to better understand the importance of their respective roles at changing operating conditions. An approach was developed which combines the simplicity and low computational and storage requests of the laminar-and-turbulent characteristic-time model with a detailed combustion chemistry model based on well-known simplified mechanisms. Assuming appropriate simplifications such as steady state or equilibrium for most of the radicals and intermediate species, the kinetics of hydrocarbons can be described by means of three overall steps. This approach was integrated in the KIVA-II code. The concept was validated and applied to a single-cylinder, heavy-duty engine. The simulation covers a wide range of operating conditions.

Combustion simulations utilizing the characteristic time combustion model have been performed for four DI diesel engines ranging in size from heavy-duty to large-bore designs. It has been found that the pre-factor to the turbulent characteristic time acts as a scaling parameter between the engines. This phenomenon is explained in terms of the non-equilibrium behavior of the turbulent time and length scales, as is encountered in the rapidly distorting, spray-induced flows of DI diesel engines. In fact, the equilibrium assumption between turbulence production and dissipation, which forms the basis for the employed k-ε-type turbulence models, does not hold in these situations. For such flows, the real turbulent dissipation time scale is locally proportional to the turbulent characteristic time scale which is determined by a typical eddy turnover time.

Intake flow simulations were carried out for a prototype DISI engine using the standard k-ε model and the RNG k-ε model. The results were compared with PTV (transient water analog) measurements. The study was focused on low load operations with engine speed at 400 rev/min. Two cases were studied, a single intake case in which one intake port was blocked and a dual intake port case. In the computations, the results show that the standard k-ε model tends to produce higher turbulence levels when turbulence is generated and decays faster when turbulence dissipates. Different turbulence models predict almost the same flow structures. However, the effects of the turbulence model on the predicted tumble and swirl ratios are significant. The TKE distributions at BDC predicted by the two models are also different. The standard k-ε model seems to be more diffusive. Good agreements with PTV data were obtained in the single valve case with the RNG k-ε model.

To model sprays from pressure-swirl atomizers, the connection between the injector and the downstream spray must be considered. A new model for pressure-swirl atomizers is presented which assumes little knowledge of the internal details of the injector, but instead uses available observations of external spray characteristics. First, a correlation for the exit velocity at the injector exit is used to define the liquid film thickness. Next, the film must be modeled as it becomes a thin, liquid sheet and breaks up, forming ligaments and droplets. A linearized instability analysis of the breakup of a viscous, liquid sheet is used as part of the spray boundary condition. The spray angle is estimated from spray photographs and patternator data. A mass averaged spray angle is calculated from the patternator data and used in some of the calculations.